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Thermal conductivity of nonlinear waves in disordered chains 
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Abstract. We present computational data on the thermal conductivity of nonlinear waves in dis- 
ordered chains. Disorder induces Anderson localization for linear waves and results in a vanishing 
conductivity. Cubic nonlinearity restores normal conductivity, but with a strongly temperature- 
dependent conductivity k(T). We find indications for an asymptotic low-temperature k ~ T 4 and 
intermediate temperature k ~ T 2 laws. These findings are in accord with theoretical studies of 
wave packet spreading, where a regime of strong chaos is found to be intermediate, followed by an 
asymptotic regime of weak chaos (EPL 91 (2010) 30001). 
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1. Introduction 

In the absence of nonlinearity (or many-body interactions in quantum systems) all eigen- 
states - normal modes (NM) - in one-dimensional random lattices with disorder are spa- 
tially localized. This is Anderson localization [1], which has been discovered fifty years 
ago in disordered crystals as a localization of electronic wavefunctions. It can be in- 
terpreted as an interference effect between multiple scattering events of the electron by 
random defects of the potential. Recent experiments on the observation of Anderson lo- 
calization were performed with light propagation in spatially random optical media 
with noninteracting Bose-Einstein condensates expanding in random optical potentials 
JUS], and with wave localization in a microwave cavity filled with randomly distributed 
scatterers Jg]. 

In many situations nonlinear terms in the wave equations (respectively, many body in- 
teraction terms in quantum systems) have to be included. Thus, a fundamental question 
which has attracted the attention of many researchers is: what happens to an initial ex- 
citation of arbitrary shape in a nonlinear disordered lattice. Nonlinearity renormalizes 
excitation frequencies, and induces interactions between NMs. Numerical studies show 
that wave packets spread subdiffusively and Anderson localization is destroyed iIt UIoIi . In 
the regime of strong nonlinearity, far from where it can be treated perturbatively, new lo- 
calization effects of selftrapping occur [11]. A theoretical explanation of the subdiffusive 
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spreading was offered in Refs. 1F7L [qL [l2Tl . It is based on the fact that the considered mod- 
els are in general nonintegrable. Therefore deterministic chaos will lead to an incoherent 
spreading. Estimates of the excitation transfer rate across the packet tail are obtained by 
calculating probabilities of mode-mode resonances inside the packet. Some predictions of 
this approach include the effect of different degrees of nonlinearity and were successfully 
tested in II 1311 . First experimental data were recently presented for repulsively interact- 
ing Bose-Einstein condensates in quasiperiodic potentials HG1 which, among others, 
confirm the basic theoretical findings. 

While the above wave packets are spreading into an 'empty' system, another set of 
related problems concerns the thermal conductivity at finite temperature T (i.e. at finite 



energy density e ~ T) Bio, 1 1711 . In the absence of nonlinearity Anderson localization 



of NMs leads to a vanishing conductivity. It was observed that even a small amount of 
anharmonicity leads to a diffusive transport of energy, and a finite thermal conductivity k 
1 1811 . Some numerical studies indicated that for low temperatures k ~ T 1 / 2 [18] which 
would indicate a singularity at zero temperatures. Other expectations claim that heat 
conductivity vanishes strictly for weak enough (but still finite) anharmonicity 11711 . 

In this report we use the theoretical results on wave packet spreading and derive a 
connection between k and T. We predict that for low temperatures k ~ T 4 which cor- 
responds to the regime of weak chaos. At higher temperatures an optional intermediate 
regime of strong chaos may lead to k ~ T 2 . We perform numerical simulations which 
indicate that our predictions are plausible. 



2. Models 

We study two different Hamiltonian models. The first model describes interacting anhar- 
monic oscillators in a quartic Klein-Gordon (KG) chain, given as 

i 

where ui and pi are respectively the generalized coordinate/momentum on the site I, and 
ti are the disordered potential strengths chosen uniformly in the interval [1/2, 3/2]. Like- 
wise, dfui = —&Hk/9ui generates the equations of motion 

3 1 . 

ui = -Qui - it, + — {ui+i + - 2ui) . (2) 

By neglecting the nonlinear terms, the KG model (Q} reduces to a linear eigenvalue 
problem. This leads to a set of NM amplitudes, A u j, with NM squared frequencies lj 2 s 
[1/2,3/2 + 4/W] and A = 1 + A/W being the width of the squared eigenfrequency 
spectrum. The NM asymptotic spatial decay is given by A v< \ ~ e~ 1 '^""' where £,(to 2 ) 
is the localization length. It is approximated [ 19] in the limit of weak disorder (W <C 1) 
as < 96W~ 2 . The NM participation number P v = 1/Xw ^-t i characterizes the 

NM spatial extent. An average measure of this extent is the localization volume V, which 
is on the order of 3.3£(0) for weak disorder and unity in the limit of strong disorder 
lEoll . The average squared frequency spacing of NMs within a localization volume is then 
d w A/V. The two squared frequency scales d and A with d < A are thus expected to 
determine the dynamics details in the presence of nonlinearity. 

Nonlinearity induces an interaction between NMs. Since all NMs are exponentially lo- 
calized in space, each of them is effectively coupled to a finite number of neighbor modes, 
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i.e. the interaction range is finite. However, the strength of this coupling is proportional 
to a characteristic energy density, e. The squared frequency shift due to the nonlinearity 
is then 5 ~ e. 

The KG model is suitable for studies of the dynamics for weak and intermediate disor- 
der W < 20, since otherwise the coupling strength between nearest neighbor oscillators 
becomes too weak and required computation times too long. In the limit of strong disor- 
der the NMs become localized on a single site, and will be coupled by nonlinearity to their 
nearest neighbors. Therefore this limit can be modeled by a chain of harmonic oscillators 
with random frequencies and purely anharmonic nearest neighbor coupling: 

nFsw = Y,Y + e i u * + \ {ui+1 ~ Ul)4 (3) 



which is an example from a class of models studied by Frohlich, Spencer and Wayne H21H 
and therefore called here the FSW model. It follows qualitatively from the KG model 
by first transforming into NM space, dropping secular terms 11211 . and finally through the 
limits of strong disorder, and of high temperatures (energies). Its quartic potential part 
however is invariant under homogeneous coordinate shifts, contrary to the KG case. 



3. Thermal conductivity 

The thermal conductivity is measured by a standard approach, when the thermal baths at- 
tached to the ends generate a temperature gradient and heat current along the chain lfl6ll . 
For the KG model we implement Langevin thermostats with two next-to-end atoms cou- 
pled to Langevin heat baths with ilx t N = —dH/dui^ + £l,jv ~ Auuv with white noise 
(£i,iv) = and (£l,jv(i)£i,iv(0)) = 2ATi i jv<$(t) where (•) denotes ensemble averages. 
T is the temperature, and A denotes the coupling strength between system and heat bath. 
For the FSW model we use Nose-Hoover thermostats adding the terms —(±ui.n to the 
respective equations of motion, where (± — u\ N /T\ n — 1. 

For the KG model the heat flux along the chain is defined as the time average of jkg = 

~2W XX^z+i + ~ u i) UM- For the FSW model the heat flux follows from the 

i 

time average of ]fsw = —\ XX^Z+l + WiX^z+i — u i) 3 UM- The thermal conductivity 

coefficient reads k = jN/ (Tjy — Ti) then lfl6ll . We make use of the mean temperature 
T = (Tat + Ti)/2 as a parameter corresponding to the energy density (e) = T, and 
(Tat - T x )/T = 0.5 further on. 



4. Predictions 

The numerical evolution of a wave packet spreading into an empty system 1F7I— TTotl yields 
subdiffusive growth of the second moment. This implies that the diffusion rate is depend- 
ing on the energy density (temperature). A theoretical approach J^HHH] postulates that 
the wave packet evolves chaotically, and estimates the diffusion rate for the KG model to 
bed 

D^ £ 2 V 2 {e) (4) 
with the resonance probability lEoll 

V{e) w 1 - e- ae/d (5) 
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Figure 1. KG chain: the size-dependent thermal conductivity kn versus N for differ- 
ent temperatures and W — 2. 



and a being an unknown constant of order one. 

For small temperatures e ~ T. Therefore, and assuming that the thermal conductivity 
is proportional to the diffusion rate ifioh . we conclude that for the KG model 

K~T 2 (l- C - bT / d Y (6) 



with b another constant of order one. In particular we find that 

k ~ T 4 , T < d (weak chaos) (7) 
k ~ T 2 , T > d (strong chaos) (8) 

Note that for larger energy densities (temperatures) the KG model enters the regime of 
selftrapping, and the above predictions do not apply anymore. The selftrapping regime 
is expected when e ~ T « 2/(3 W) fl 1 Oil . It is the asymptotic high temperature limit of 
the KG model, where the thermal conductivity should drop with increasing temperature 



01711 . Thus, depending on the strength of disorder, we expect that the thermal conductivity 
grows with temperature as T 4 , then (for not too strong disorder) it grows as T 2 , and finally 
goes through a maximum and starts to decrease with further increase of temperature. The 
three phases are summarized in the inset of Fig. 3 in Ref. II 1 CHI . 

The application of the above approach to the FSW model leads to d ~ 1 and to the 
absence of a selftrapped regime. Instead, for large temperatures, the leading order an- 
harmonic potential terms of the FSW model bring it close to the momentum conserving 
/3-Fermi-Pasta-Ulam chain (FPU). The FPU chain has a thermal conductivity which di- 
verges with increasing system size iflfjll . Therefore k will increase with system size also 
for the FSW model for large temperatures, until the harmonic oscillator potential parts 
can not be neglected anymore, leading to a final saturation of the thermal conductivity 
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with further system size increase. Indeed, the harmonic oscillator potential parts are vi- 
olating total momentum conservation, and therefore the thermal conductivity stays finite 
for large enough system size H 1 OL 1 1 '711 . Since the momentum conserving terms of the FSW 
model become more dominant at larger temperatures, it follows that the higher the tem- 
perature, the larger the system sizes are needed, and the larger the saturated conductivity 
value. The crossover to this pseudo-FPU regime can be expected when the temperature is 
of order unity (since all relevant FSW model parameters are of order unity). Therefore we 
conclude, that we will not see a pronounced regime of strong chaos, instead the regime of 
weak chaos at low temperatures should gradually transform into the pseudo-FPU regime. 

5. Numerical results 

The thermal conductivity is computed along the above lines, with additional averaging 
over disorder realizations. 

5.1 KG chain 

We first present results for the disorder strength W — 2. In that case the localization 
length £ m 25 and the localization volume V « 75. The thermal conductivity kn is 
plotted for various temperatures versus system size N in FigQ] For large temperatures 
T > 0.1 the conductivity reaches its saturated level n at system size N ks £. For the linear 
wave equation, exponentially localized NMs induce a dependence nu near ~ e~ N ^. In 
the presence of nonlinearity the exponential decay with N will stop at the saturated value 
of k. Assuming that n ~ T a it follows that the system size for saturated conductivity 
values scales as N ~ — a£ In T. Thus for lower temperatures the system size needed to 
reach saturated values of n increases. This is due to the dropping of the saturated value 
of k with T. We were able to reach temperatures as low as T = 0.005 and conductivity 
values k = 0.008. According to Ref. II 1011 (inset in Fig. 3) the regime of strong chaos is 
realized for 0.01 < T < 0.3. Therefore we expect to observe « ~ T 2 for 0.01 < T < 1, 
followed by a maximum and a further dropping with increasing temperatures. These 
predictions are indeed compatible with the numerical results in Fig|2] We also compute 
the conductivity for a reference ordered chain with £/ = 1, which coincides with the high 
temperature part of k(T) for the disordered case. Therefore the observed maximum is 
due to disorder, as well as the decrease of k with further lowering of the temperature. 
The low temperature data are also consistent with the T 2 law. However the temperature 
window is too narrow to conclude about a confirmation. Also we expect that at even lower 
temperatures (so far not reachable with our computational means) the conductivity should 
cross over into the weak chaos regime and follow the T 4 law. 

In order to find evidence for the asymptotic low-temperature behavior k ~ T 4 we chose 
a stronger disorder W — 6. In that case £ ss 3 and V « 10. The conductivity is reduced 
by two orders of magnitude. We are therefore quite limited in the temperature window, 
and do not obtain conclusive data (see Figf5]l. We again find a maximum in the k(T) 
dependence, and a clear indication that the small temperature conductivity drops faster 
than T 2 with temperature. Yet we do not have enough low temperature data to conclude 
about the exponent. 

5.2 FSW chain 

To observe the asymptotic low temperature behavior of the conductivity, we explore the 
strong disorder limit and study the FSW chain. Since the localization length £ = for 
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the FSW chain, we can use small system sizes for small temperatures, while (as discussed 
above) we will need larger system sizes to obtain the saturated conductivity values for 
larger temperatures. In Figf3]we plot the conductivity n versus temperature for different 
system sizes N = 16, 64, 256. Indeed, we observe a significant difference between the 
data for temperatures T > 1, while the data agree well for lower temperatures. Therefore 
we can reach lower values of the conductivity as compared to the KG chain case. Yet for 
k < 10~ 4 (which corresponds to T ~ 0.01) we again reach our computational limitations, 
since the number of realizations and the integration times, which are needed in order to 
reduce the impact of fluctuations, increase drastically. The numerical results in Figf3] 
indicate that for low temperatures the conductivity behaves as k ~ T 4 . Also we observe 
at larger temperatures a crossover to the predicted pseudo-FPU regime. Again we need 
more data in the low temperature regime to be more conclusive about the conductivity 
data. 



6. Discussion and outlook 

We computed the thermal conductivity for nonlinear disordered chains and obtained re- 
sults which indicate that theoretical predictions based on wave packet spreading are plau- 
sible. In particular we find indications for an asymptotic low-temperature n ~ T 4 and 
intermediate temperature n ~ T 2 laws. To be more conclusive, we need data at even 
lower temperatures than those presented here. The needed CPU times are however cur- 
rently exceeding our capabilities. Possible solutions include the usage of massive parallel 
computing equipment, or the search for other models which will allow to make the needed 
observations with standard computational equipment. 
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Figure 3. FSW chain: k(T) for system sizes N — 16 (filled circles), N = 64 (filled 
squares), and N = 256 (filled diamonds). The dashed line corresponds to the law T 4 . 
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